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Abstract 

We discuss the complex phase generated in three pion correlation functions. The lowest order contribution 
to the phase is of order q 2 R/K, where q is a typical relative momentum, K is a typical center of mass 
momentum and R is a typical radius parameter. This contribution is of purely kinematic origin. At next 
order we find a generic contribution of order (qR) 3 which is a result of odd modifications to the source emission 
function. We argue, that the scale for typical HBT correlations in ultrarelativistic heavy ion collisions is 
q/K <C qR ~ 1, so that the third order correction actually dominates the phase in the experimentally relevant 
momentum range. We study in detail such contributions which arise from source asymmetries generated by 
flow, the source geometry and resonance decays. 



I. INTRODUCTION 



Recently the first study of three pion correlations (7r + 7r + 7r + ) in heavy ion collisions was reported by 
NA44 Q. The quality of the data was rather limited due to statistics, but we can expect a steady increase 
in statistics in forthcoming experiments. This should allow us to explore large parts of the phase space of 
the three pion correlation function in the foreseeable future. 

The theoretical understanding of the three-pion correlation function is still rather limited. In an earlier 
paper |2| we constructed the three pion correlation function under the assumption that the source was inco- 
herent and only interfered through Bose-Einstein correlations. Thus the three-particle correlation function 
could be derived from the two-body correlation function except for a complex phase, which will be studied 
here. Any deviation found experimentally from this prediction would have signalled new physics. We found 
satisfactory agreement with the available data, but due to the poor statistics could not make any conclusive 
statements. 

Heinz and Zhang || explored the structure of the complex phase using an expansion in terms of the 
relative momentum of the two emitted particles. They found the lowest order contribution to be in the order 
q 2 R/K. Here, the relative momentum is given by q = (ki — k2) the average momentum by K = (ki + k2)/2 
and ki is the single particle momentum of particle i. This contribution is of purely kinematic origin. Typical 
heavy ion sources in nuclear collisions are of size R ~ 5 fm, so that interference occurs predominantly when 
q^h/R ~ 40 MeV/c. Since typical particle momenta are fc, ~ K ~ 300 MeV, we find that q/K -C 1 and the 
lowest order contribution to the phase is very small. 

In this paper we will extend the two calculations mentioned above and explicitly investigate the modifi- 
cations due to generic three pion correlations. In section 2 we review the general derivation of the phase in 
3 particle correlations in terms of an expansion up to cubic terms in the relative momenta q. In section 3 
we estimate the size of the phase for a number of physical effects known to be present in relativistic heavy 
ion collisions and which produce asymmetric sources leading to a complex phase in the 3 particle correlation 
function. These are flow, resonance decay and asymmetric source geometry. We find some contributions 
to be sizable, but the overall effect off such a phase on the measurable correlation function turns out to be 
rather small and any detection will require exceedingly good resolution in the regime for which qR ~ h. 
Finally, in the conclusions, we summarise our results and discuss some experimental consequences. 



II. GENERAL FORMALISM 

The three pion correlation function for incoherent sources is given by [§,f§-|§ 

C 3 (ki, k 2 , ks) = 1 + F\2 F21 + F 23 F32 + F31 Fi3 + F12 i*23 F3I + F21 F 32 F13 

= 1 + \F 12 \ 2 + \F 23 \ 2 + \F 31 \ 2 + 2 Re [F 12 F 23 F 31 ] , (1) 
where Fij is the Fourier transform of the source emission function S(x, K) [Q 

Fij ee F(q id , Kij ) = ffoSfrKtfaxpdqyz) = ^ (JJ) 



yfj d 4 x S(x, h) J d 4 x S{x, kj) 



We used here the relative momentum — ki — kj and the center of momentum variable Kij = (ki + kj)/2. 
The emission function S(x, k) is the probability of emission for a pion from space-time point x with momen- 
tum k. It is related to the experimentally measured pion single particle spectrum 

E,^t=[ d A x Six, k) , (3) 



where Ek = fco = Vk 2 + m 2 is the on mass-shell energy. 
The relative momenta of three particles satisfy the relation 



<7l2 + <?23 + 931 = , (4) 

i.e., they span a triangle. This automatically assures translational invariance of the 3-body correlation 
function. Any translation in space-time by a distance xq will lead to an extra phase factor exp{iqij ■ xq) in 
Fij but Eq. (H) insures that the triple product of the three phase factors cancel in the 3-body correlation 
function ([!]). 

The source emission function in the numerator of Eq. (Q) is not evaluated for a momentum fej , but rather 
for the center of mass momentum . Since all particles are detected on shell, we will have to evaluate the 
center of mass momenta in (|^) slightly off shell 

K « = E « ( 1 + s4 + °(^)) ' (5) 

This on-shell constraint for the two detected particles pushes the source emission function slightly off-shell 
by an amount 

S(x,K° j ,K ij )=S(x,E K ,K ij ) +J^^L-( x ,E K ,K ij )+0(t$ j ) . (6) 

This correction is of order (q/K) 2 and is also present in the 2-particle correlation function. Since the off-shell 
structure of the source emission function is not accessible to us we will neglect this contribution, like in the 
2-particle case, with the remark that a possible signal could be due to this correction. 

With this in mind, we can rewrite Eq. (Q) using on shell variables Fij — -F(qij, K^) and explicitly depict 
the phase 0y = <£(q i: ,-,Ky) 

C 3 (k 1 ,k 2 ,k3) = l + |F 12 | 2 + | J F 1 2 3| 2 + |F 31 | 2 + 2| J F 1 12 | \F 23 \ \F 31 \ cos (0 12 + 23 + fa) . (7) 

All information about possible imaginary contributions to F are now contained in the cosine of the 3 phases, 
4>ij, defined as 



tan<fe = r-^ ■ (8) 

' Re iV 



The phases arise from the triple product of Fourier transforms in Eq. For the following discussion it 
is important to realize, that the three particle correlations as given above are defined over a 9 dimensional 
momentum space. This space can be either described by the 3 vectors hi , k 2 , k 3 , or preferably by the center 
of momentum of the three emitted particles K = (ki + k 2 + k 3 )/3 and two relative momenta, like qi 2 and 
q 23 . The kinematic transformation into center of momentum and relative momenta will actually fix the 
lowest order contribution in our expansion below. 

To determine the phases we have to extract both the real and imaginary pieces of the Fourier transforms 
of the source emission function. The source emission function itself is real. The only source for an imaginary 
contribution is thus the Fourier transform via the exponential of the relative momentum. We introduce the 
symmetric and antisymmetric part, S s and S a , of the source emission function by 

S s (x - (x), K) = i [S (x - (x) , K) + S (-(x - (x)) , K)] 

S a (x-(x),K) =^[S(x-(x),K)-S(-(x-(x)),K)} . (9) 
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The average of an operator £ is defined as 



(i) = 



Jd 4 xiS{x,K) 



(10) 



Jd i xS(x,K) ' 

The space-time integral over the asymmetric part is identical zero, while the space-time integral over the 
symmetric part provides the normalisation of the source emission function. 
The phase <f>ij in Eq. (||) is now 

tan ^ _ 1 J (fix S a (x,E K ,Kij) expjiqjjx) 
13 i J d 4 x S s (x, E K ,Kij) exp(iqijx) 

which clearly demonstrates that the odd space-time moments of the source emission function S generates 
the phase. If S is even, then the Fourier transform in (^|) will be real and the phases 4>ij vanish. 
If we expand the exponential in Eq. (11) for small values of qij 



exp (iqijx) = 1 + iq tj x - - (qijx) 



{qijxf 



0{q%) 



we obtain for the phase 



tan<fe = (fax)) ~ - {{ qij {x-{x))f)) + 0{q%) 



(12) 



(13) 



The sum of the 3 phases is then 



1 



^13 



»23 - 2 9l2<?23 



d{x,j) d{x v ) 



r M v A 1 U v A 1 
^ L9l2?12923 + 923923912J 



d 2 (x, 



d 2 (x v ) 



d 2 (x x ) 



dK v dK x 8K x dK^ dK^dK" 



1 



ttfe + 92 3 ) A ( (X - (!))„ (I ~ (X)) v (X - (x)) x ) + 0(q%) 



(14) 



This result was also obtained by Heinz and Zhang ||. To obtain this formula we made use of the triangle 
relation (^J) which assures that the three linear terms sum up to zero. The on mass shell constraint in (^) 
fixes the time components of the four vectors. These redundant components are not explicitly eliminated. 

The last term in Eq. (|l4|) is of order (qR) 3 and is generic. In contrast the first two terms are a result 
of choosing the three momenta K, q\2 and 523 to span our 9 dimensional coordinate space and are of order 
q 2 R/K and q 3 R/K 2 respectively. Since typical particle momenta are fcj ~ K ~ 300 MeV and typical heavy 
ion sources in nuclear collisions have a size R ~ 5 fm, we find that q/K <C qR ~ fi, so that the generic 
(<7-R) 3 -contribution can actually dominate in the experimentally relevant momentum regime. 

In a simple but realistic model for flow, we will demonstrate in the next section why contributions of order 
q 2 R/K are so small. Afterwards we check in two further models for resonance and source geometry if we 
actually can produce a significant phase effect at scales of order qR ~ h. 



III. ASYMMETRIC SOURCES 

From Eq. ( |l3| ) and (|l4|) we see that there are contributions to the phase due to odd space-time moments 
of the source emission function of order q 2 R/K and (qR) 3 . To evaluate these contributions we need a model 
for the source. In the following we will study a number physical effects, which are known to be present in 
relativistic heavy ion collisions, and which result in asymmetric sources. The examples we will study in the 
next four subsections are flow, resonances, moving and bursting sources respectively. We will estimate their 
quantitative influence on the phases ^ij ■ 
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A. Flow 



The most common ansatz for particle production and collision dynamics in ultrarelativistic heavy ion 
collisions is the Bjorken scenario. One assumes cylindrical symmetry and requires local thermal equilibrium 
with longitudinal Bjorken flow (u z = z/t) as well as transverse flow v through a Boltzmann factor. Thus 

S(x,K)~e- K - u / T S x (x). (15) 

where u = 7(u)(v, sinh^), cosh(?y)) is the flow four-vector so that 

K ■ u — m^7(t>)(cosh(7/ — Y) — (3k ■ v) . (16) 

Here, r = \Jt 2 — z 2 is the invariant time, r\ = 0.5 ln(i + z)/(i — z) the space-time rapidity, and v = (v x , %) 
the transverse flow. The pair velocity is /3k = K/A°. The transverse flow contribution in Eqs. (^) and 
( pf ) provide a strong ^-dependence. Further if -dependences in the emitting source, S x (x), would intro- 
duce additional A-dependences. All these A-dependences should contribute to the second order kinematic 
correction in Eq. (H). 

When transverse flow is present only a few terms contribute. They all vanish though, if we boost into 
the so called longitudinal center of momentum frame (LCMS). In this frame we have Y — and the center 
of momentum of the pair is parallel to the x- direction. This direction is labeled the outward direction (o), 
while the beam axis or z-direction is labeled the longitudinal direction (1). Perpendicular to these is the 



sidewards (s) or y-direction. With this in mind we can evaluate the (q/K) 2 contribution of Eq. (14) 



d{x ll ) d{x v ) 



= <7i2<?23 \{x^u u ) - {x v uj\ = . (17) 



In the last step we made use of the fact, that (y) = due to cylindrical symmetry and reflection symmetry in 
the x — z plane. We also boosted to the LCMS in which Y — 0, so that (z) = (u z ) = and v = (v x , 0) points 
into the outward direction. Furthermore, due to symmetry any choice for v x has to be an even function in y 
and z, so that (v x z) = (v x y) — 0. Finally, in the LCMS the pair velocity simplifies, so that = f3 x qfj- As 
a result the xt— and tx— contribution in Eq. ( |l7j ) cancel each other. 

This result demonstrates, that the standard A-dependent contribution to the source emission function, 
i.e., flow does not produce sizable contributions to the kinematic correction of second order. This correc- 
tion actually vanishes in the LCMS. We will see in the following examples that higher order, but generic, 
contributions like resonance and source geometry can in principal dominate the experimentally accessible 
momentum range. 

B. Resonances 

We can study the influence of resonances on the phase within a simple model S . Under the assumption 
of classical propagation, the resonance travels an extra distance Ax = u r At before it decays and produces a 
pion. At is hereby the life time of the resonance and u r its velocity. If, furthermore, the resonance life time 
is exponentially distributed with a decay width T r = l/r r we obtain after averaging over the resonance life 
time a modified source emission function. Its Fourier transform reads 

F r (q.ij, Kij) - J d A x S(x, E K , K y ) exp^^x) (1 - iq^ ■ w r r r ) _1 . (18) 
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Every resonance r will supply such a contribution to the correlation function and x refers to the space-time 
production point of the resonance. 

We can evaluate the phase in (||) for a source consisting of one resonance only. For a resonance velocity 
independent of the space-time production point x we find 

tan^jj = u r ■ qij r r . (19) 

We would like to investigate the significance of this result and compare it with current experiments. Due 
to the limited statistics in the experiments one reduces the 9 dimensional momentum space, on which the 
cosine defined in Eq. (Q) depends, down to one invariant momentum 

Ql = <7l2 + 931 + 923 . (20) 

which is used to analyse the data. 

In a previous paper we evaluated the radial part of the Fourier transform F, neglecting any phase contri- 
bution and found that the correlation function in Eq. (0) can be rewritten as 

2 

C 3 (Q 3 ) = 1 + 3 A 2 exp(-y) (1 + 0(x 4 )) 

x 2 

+ 2 A*' 5 exp(-y ) (1 + 0(x 4 )) cos (0 i2 + 23 + hi) , (21) 

where x = Q 3 R av and R av is an average source emission size parameter. R av is obtained by forming the 
average mean square of the experimentally measured radii from two particle correlations. The parameter A 2 
accounts phenomenologically for a number of effects like a partially coherent sources, long lived resonances, 
final state interactions and Coulomb screening effects. These effects tend to reduce the value of this parameter 
form its ideal value 1. It's value is also taken from the experimentally determined two particle correlation 
functions. The corrections of order x A are due to the non-spheriocity of the source emission function. These 
corrections are found to be rather small experimentally, in the order of a percent, i.e. the source is close to 
spherical. 

In calculating the phase we consider for simplicity only the momentum component in detector or In- 
direction, i.e. the outward component Qij.o, and set all other momentum contributions to zero. With this 
simplifications we can evaluate the cosine in terms of only the outwards momenta 



(COS (012 + 023 + 03l))(Q3,o) = 

/ dq 12 ,o <f?31,o dq 23 . COS (0i 2 + 023 + 031 ) HQ 3 ,o ~ 9l2,o ~ ?31,o ^ 923,o) %12,o + 931,0 + <723,o) 



(22) 



/ dq 12 , dq 31:0 dq 23yQ S(Qj - q\ 2 - q\ lo - g| 3 o ) 5{q X2 ,o + <?3i,o + 923,o) 

The ^-function of the sum of the relative outward momenta assures that the 3 particles span a triangle. 

In figure 1 we plot the cosine defined in Eq. ( |22] ) for different resonances. At large relative momentum 
each phase approach ir/2 as seen from Eq. (|l9|). Therefore the cosine of the three phases vanish for large Q 3 , 
on a scale set by the decay width l/r r . In the momentum region where correlation functions are sizeable, 
Qij^bO MeV, the short lived resonances as the K*, A and p do not produce a sizeable phase. The long lived 
resonances as 77, rj' and Kg cannot be seen at all due to their small form factor ~ (1 + u r ■ qr r )~ 2 which 
effectively removes these long lived resonances from Bose-Einstein correlations and leads to an effectively 
smaller A Only the lo could have a significant influence on the phase at momenta Qy ~ 50 MeV. However, 
it would be small due to its own form factor and due to the small fraction of pions that are decay products 
of w's. 
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FIG. 1. The phase contribution for different resonances as a function of the momentum Ql j0 = Qi2, + 913,0 + 923,0 • 
The fall-off of the curves scales with the resonance life time. The transverse momentum of the pion pair is chosen 
#i = 0.7. 



The important part to notice is, that the phase contribution is suppressed additionally by a Gaussian form 
factor as in Eq. (|2l|). Even if the cosine is sizeable, like for example for the to, the overall modification to the 
measured correlation function is still rather small due to the form factor suppression. Only an experiment 
with excellent resolution will be able to resolve the resonance contribution. 



C. Moving surfaces 

To have a detectable contribution to the phase we need to create an odd modification of the source emission 
function, which appears at a momentum scale small enough, so that it is not suppressed to drastically by the 
Gaussian form factor. On the other hand such a contribution should be weighted strong enough to overcome 



the suppression for small momenta shown in Eq. (13). In this section we try to construct an odd source 
geometry, which has this properties. 

Inspired by hydrodynamical models we assume a source which predominantly emits particles from a thin 
surface layer. Such a source emission function can be described as 

S A (x)~6(R(T)-r±), (23) 



where r± is the transverse radius and r is the proper time. The source is cylindrically symmetric with a 



radius R(t) = Rq (1 — (t/t/) q ). The transverse radius starts out at i?o and the source disintegrates at time 
T f- 

It is straightforward to evaluate the phase such an extreme geometry produces. We find 

. (- 1 )" q2 " +3 

a cos a - sin a + 2^„=o ( 2 n+i)! (2n+ Q +3) 
tan(0 a ) = : (-d^ +s- , (24) 

COSO + O SUIO - 1 - £„ =0 (2»)! (2n+a+2) 

if we assume, that the outward component of the momentum is dominant, i.e. if we neglect the other two 
components. The only remaining variable is a = (3j_qij^ Tf, where (3± is the transverse velocity of the pion 
pair. For a = lwe find </>i = —a/2 and in the limit of a — > wc obtain 

0a=~|«|^+O(« 3 ). (25) 

6 6 + a 

The term linear in a vanishes once the cosine of the sum of the three phases is evaluated due to the triangle 
constraint, Eq. (^). We plot in figure 2 the difference A</> between the phase given in Eq. ( ]24| ) and the linear 
term of Eq. (p5|). This difference is a direct measure for the strength of the cubic correction in Eq. (O). 




FIG. 2. The difference A<j> between the phase 4>a and its limit for small q in the case of surface emission. 

The contribution for a — 1 to the phase difference vanishes identically. A surface, moving with constant 
speed does not produce any phase. In this case the product tR(t) oc t (t/ — t) is symmetric around 
Tf/2. The additional factor of r in the product comes from the integration measure suitable for the Bjorken 
scenario. For accelerated surfaces, where a/lwe find a sizable decrease at values fi±Tfqo ^> 5. For typical 
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values of the transverse momentum 0± ~ 0.7 and typical time scales t/ of a few Fermi this corresponds to 
outward momenta in the GeV region. We have to rule out such a surface process. 



D. Bursts 

Another possibility to geometrically obtain odd modification functions to the source emission is asymmetric 
bursts 

SvAx) ~\ [(1 - e) S(t - t c ) + (1 + e) S(t -t c - At)} S s (r). (26) 

This source corresponds to two bursts of particles emitted at times t c and t c + At from a spatial emission 
source function, S s (r); the latter is arbitrary and irrelevant as long as it is symmetric. The two contributions 
are weighted by e such that |e| < 1. 

The phase for such a source is from Eq. (|§) 

tan - qi(t c + At/2)) = e tan (<j 4 At/2) . (27) 

where = Ei — Ej ~ [3±q - The contribution proportional to the mean time q^ (i c + At/2) on the right hand 
side cancels when the three phases are added due to the triangle constraint (^). Thus the phase effectively 
vanishes for equal weights e = 0, if one weight vanishes e = ±1 or for zero time separation of the two flashes. 
This is simply because the source is symmetric in all these cases. The scale of variation of the phase is given 
by the At and the size of the phase is given by the asymmetry parameter e. 




FIG. 3. The phase contribution in the case of volume emission from two instantaneous sources separated by a 
time interval At. The weight e is set to .25, .5 and .75. 

In figure 3 we plot the cosine of the three phases according to Eq. ( [22] ) for different values of e, assuming 
again that the outward component of the relative momentum is dominant. The phase scales in the variable 
/3j_Q3, Ai. For outward momenta close to (3±Qs t0 At/2 ~ 2 one of the phases in Eq. ( p7| ) changes sign and we 
see a strong signal. In the experimentally accessible momentum range this would correspond to a temporal 
separation of the bursts in the order of 10 fm/c. While such a scenario is rather unrealistic it might open in 
the long run some new ideas and approaches to investigating possible phase signals. 

IV. CONCLUSION 

We have studied 3-body correlations for incoherent sources which can be calculated from 2-body corre- 
lations except for a phase. The phase is due to odd space-time moments of the source emission function 
and vanish at small momentum transfer to order (qR) 3 and q 2 R/K, where the former dominates in the 
experimentally relevant momentum regime in relativistic heavy ion collisions. 

Effects of flow, resonances, moving and bursting sources were studied as they are known to be present in 
relativistic heavy ion collisions, and result in asymmetric sources. Their influence on the phases <pij was only 
significant at large relative momenta where the form factors, F(q), and thus also correlation functions were 
small or in some extreme scenarios, where all particles where emitted in 2 sudden bursts. In both cases, the 
experimental extraction of the phases needs very high resolutions. The phase in 3-body correlation functions 
seems to be very elusive. 
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FIGURE CAPTIONS 

Figure 1. The phase contribution for different resonances as a function of the momentum Q\ a = 
Q12 o + I13 o + 923 o- The fall-off of the curves scales with the resonance life time. The transverse mo- 
mentum of the pion pair is chosen (3±_ = 0.7. 

Figure 2. The difference A<j> between the phase 4> a and its limit for small q a in the case of surface emission. 

Figure 3. The phase contribution in the case of volume emission from two instantaneous sources separated 
by a time interval A. The weight e is set to .25, .5 and .75. 
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